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ABSTRACT 

It has been shown that HD molecules can form efficiently in metal-free gas collapsing 
into massive protogalactic halos at high redshift. The resulting radiative cooling by HD 
can lower the gas temperature to that of the cosmic microwave background, Icmb = 
2.7(1 + z)K, significantly below the temperature of a few xlOOK achievable via H2- 
cooling alone, and thus reduce the masses of the first generation of stars. Here we 
consider the suppression of HD-cooling by UV irradiation in the Lyman- Werner (LW) 
bands. We include photo-dissociation of both H2 and HD, and explicitly compute the 
self-shielding and shielding of both molecules by neutral hydrogen, HI, as well as the 
shielding of HD by H2. We use a simplified dynamical collapse model, and follow the 
chemical and thermal evolution of the gas, in the presence of a UV background. We 
find that a LW flux of J cl -it,HD ~ 10~ 22 erg cm~ 2 sr -1 s _1 Hz -1 is able to suppress 
HD cooling and thus prevent collapsing primordial gas from reaching temperatures 
below ~ 100K. The main reason for the lack of HD cooling for J > J cr it.HD is the 
partial photo-dissociation of H2 , which prevents the gas from reaching sufficiently low 
temperatures (T < 150K) for HD to become the dominant coolant; direct HD photo- 
dissociation is unimportant except for a narrow range of fluxes and column densities. 
Since the prevention of HD-cooling requires only partial H2 photo-dissociation, the 
critical flux J cr it,HD is modest, and is below the UV background required to reionizc 
the universe at z ~ 10 — 20. We conclude that HD-cooling can reduce the masses of 
typical stars only in rare halos forming well before the epoch of reionization. 
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The first generation of stars are believed to be much more 
massive (~ IOOM0 ) than typical stars i n stellar populations 
in the low-reds hift universe (~ IMq ; iBromm et al.l |2002| ; 
I Abel et"aL 2002). This has many important consequences in 
the early universe, for reionization, metal-enrichment, the 
formation of seed black holes at very early times, and the 
observability of first-generation galaxies. 

The high masses result from the thermodynamical prop- 
erties of H2 , the main coolant in low-temperature gas with a 
primordial composition. In particular, H2-cooling becomes 
ineffective at temperatures below ~ 200K. HD molecules, 
can, in principle, cool the gas to much lower temperatures, 
but until recently, the abundance of HD in the early universe 
was believed to be too low for it to be important. 

It has recently been pointed out that significant HD 



can form in metal-free gas, due to non-equilibrium chem- 
istry, provided that the gas has a large initial electron frac- 
tion. This can occur, for example, in "fossil" gas that was 
ionized by a short-lived massive star, prior to it being ex- 
tinguished, or in collisionally-ionized halos with virial tem- 
peratures above « 10 4 K. It has been shown that the result- 
ing radiative cooling by HD can then lower the gas tem- 
perature to values near that of the cosmic microwave back- 
ground, T C mb = 2.7(1 + z)K, i.e. to ~ 30K atz~ 10. This 
would decrease the expected masses of the stars that form 
in ionized halos by a factor of ~ 10 below that which is 
possi ble if HD-cooling is neglected (e.g. Ijohnson fc Bromml 
2006). Q Thus, a second mode of star formation has been 
proposed, giving rise to Pop. III. 2 stars that can form 
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1 Using three-dimensional simulations, iMcGreer fc Brvanl feOOSft 
found that HD lowers the expected Pop. Ill masses less dramati- 
cally, but still has an important effect. 

2 Adopting the terminology suggested by IBromm et alj (2009 ) . 
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as soon as a small number of Pop. III. 1 stars have ini- 
tiated the epoch of reionization, and whose masses are 
only a few tens of s o lar m a sses (TNakamura fc Umemural 
20021; iMachida et alj 120051; iNagakura fc Omukail 120051; 



Johnson fc Bromm 20061. iRipamontT 2007 ; Yoshida et alj 
2007l a: lYoshida et alJl2007l b). 



These conclusions could potentially be revised, however, 
due to the suppression of HD-cooling by UV irradiation of 
the gas cloud. Alth ough this possibility has been raised in 
the l iterature (e.g. I Johnson fc Bromml 120061 ; lYoshida et al.l 
l2007l b). previous work has not included a detailed treat- 
ment of of the impact of UV irradiation on HD-cooling, 
including photo-dissociation of HD by radiation in its Ly- 
man and Werner (hereafter LW) bands, taking into account 
the shielding that occurs in the optically thick regimes. Such 
UV radiation will exist in the early universe, and can sup- 
press H2-cooling in low-mass halos at high redshifts (e.g., 
lHaiman et al.l ll997). The main goal of this paper is to as- 
sess whether HD-cooling can be similarly suppressed by UV 
radiation, and to compute the critical UV flux for the HD- 
destruction. 

In order to do this, we perform "one zone" calcula- 
tions with a simplified density evolution, while following 
the gas-phase chemistry and thermal evolution of the gas, 
including the impact of H2- and HD-dissociating LW ra- 
diation. In general, collapsing gas clouds become optically 
thick to this radiation, so that the effects of self-shielding 
are non-negligible. Our treatment includes self-shielding of 
HD and H2, shielding of both species by neutral hydrogen 
(HI), and shielding of HD by H2. We provide useful fitting 
formulae for these shielding factors, analogous t o the case of 
H2 self-shielding studied bv lDraine fc Bertoldil l| 19961 ) (here- 
after DB96). 

The rest of this paper is organized as follows. In §[2] we 
describe our chemical, thermal, and dynamical modeling. §[3] 
presents our results on the critical flux required to suppress 
HD-cooling, followed by a brief discussion of the potential 
cosmological implications and primary uncertainties in § [4] 
We summarize our main results and offer our conclusions in 
§ [5] Throughout this paper, we adopt the standard ACDM 
cosmological background model, with the following pa- 
ramet ers: Rdm = 0-233, O b = 0.0462, fi A = 0.721, and h = 
0.701 |Komatsu et alj|2009l) . 



comparable that of the intergalactic medium (IGM) at high 
redshift, n ~ 10~ 7 (1 + z) 3 cm -3 , and temperature T sa 10 4 
K. We model the second case by a pre-imposed density evo- 
lution obtained from the spherical collapse model. 



2.1 One-Zone Spherical Collapse Model 

We adopt the model for homologous spherica l collapse that 
has b een used in several previous studies (e.g. lOmukai et ail 
2008; hereafter OSH08). This simple one-zone treatment 
prescribes the density evolution of the baryonic and dark 
matter (DM) components of a collapsing halo. Both are ini- 
tialized with zero velocity at the turnaround redshift, set 
throughout this paper to z = 17. The density of the in- 
falling gas evolves on the free-fall timescale and that of the 
DM is given by a top-hat overdensity until virialization, after 
which it remains constant at its virial value. Compressional 
heating is included in the thermal model, along with the 
processes listed in § [2j2. 

Unless stated otherwise, we take the radius of the cloud 
to be R c = Aj/2, where the Jeans length is given by 



Aj = 



7rfcsT gas 

G/9gas/im p ' 



(4) 



Here ks is Boltzmann's constant, T gas is the gas tempera- 
ture, /J, is the mean molecular weight, and m p is the mass 
of the proton. Note that the size of the cloud is required, in 
practice, only in our calculations of the self-shielding factors 
(see below), in order to specify the column densities of H2, 
HD, and HI. 

While this model is a vast simplification of the physics 
of a collapsing halo, it nonetheless has been shown to 
mimic the thermal and chemical evolution seen in full three- 
dimensional hydro dynamical simulations very well (see, e.g., 
IShang et al 1 l2010l - hereafter SBH10 - for a direct compar- 
ison). The exception is the shock-heating that occurs in the 
early stages of collapse and is not present in the one-zone 
model, which prescribes a smooth "free-fall" evolution. For 
a detailed description of the spherical collapse model, the 
reader is referred to the recent work by OSH08 and refer- 
ences therein. 



2 MODEL DESCRIPTION 

The formation of HD o ccurs primarily thro ugh the following 
reaction sequence (e.g. iGalli fc Pallall2002i ): 

R + e~^H~ + hv (1) 
H + H"^H 2 +e" (2) 

D+ + H 2 ^HD + H + . (3) 

Thus, in order to form a significant abundance of HD, a large 
initial electron fra ction is required to catal yze the forma- 
tion of H2 (see, e.g. Ijohnson fc Brommll2006l . and references 
therein). In primordial gas this can be achieved by photoion- 
ization (e.g. by short-lived Pop III. 1 stars), or by collisional 
ionization (in sufficiently massive halos). We model the first 
case by a constant, low-density gas, initially at a density 



2.2 Chemical and Thermal Model 

We model a gas of primordial composition using a re- 
action network which comprises 47 gas-phase reactions 
amongst the following 14 chemical species (and photons): 
H, H + , H~, He, He+, He 2+ , H 2 , Hj, D, D+, D", HD, 
HD + , and electrons. Our choices for the selection of species 
and their initial a bundances are conventional (see, e.g., 
IGalli fc Pallal [ 19981 ). but we do not include any lithium 
species or other potential coolants (e.g. H ^~), as they con- 
tribu te very little to the total cooling (e.g. iGlover fc Savinl 
2009) and are not important in the context of this paper. 

2.2.1 Hydrogen and Helium Chemistry 

The collisional rate coefficients for reactions among hydro- 
gen and helium species only, and cross-sections for photo- 
ionization, are taken from the recent compilation by SBH10. 
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However, the rate for H2 photo-dissociation (&28 hi the afore- 
mentioned compilation) is modified to fcdiss,H 2 = 1-39 x 
10 -12 x /3 x /shield in order to match the optically thin rate 
we calculate (see § I2.3|1 . The total shielding is parameterized 
by a shield factor, /shield (see below), and the rate is normal- 
ized by the parameter /?, as described by in Appendix A of 
OM01, which specifies the intensity of blackbody radiation 
at the average LW band energy (12.4 eV) relative to that 
at the Lyman limit (13.6 eV). For the two spectral types 
we consider (described in § |3J, /3 = 3 for the T4-type and 
P = 0.9 for the T5-type. 



2.2.2 Deuterium Chemistry 

The chemical network includes 19 reactions involving the 
five deuterium species, for which w e use the collisional rate 
coeffic ients from the compilation bv lNakamura fc Umemural 
(2002). However, we replace the rates D2, D3, D7, and D9 
given therein (as referenced i n the s ource) with the corre- 
sponding updated rates from ISavi n| (120021) (for the charge 
exchange reactions, D2 and D3) and Galli fc Pallal (|2002T ) 
(for D7 and D9). The HD photo-dissociation rate is given 
in § 12.31 and is normalized with the /3 parameter in the same 



manner as described above for the H2 photo-dissociation 
rate. 

We take the cosmological D/H ratio to be 4x10 
by number, following recent studies on HD -cooling 
(e.g. I Johnson fc Bromml 120061, lYoshida et al.l|2007h and in- 
spired by the model of iGalli fc Pallal 1)19981 ) , which pro- 
vides a value of D/H = 4.3 x 10 -5 . This adopted value 
is likely overgenerous for the primordial deuterium abun- 
dance, however, in light of recent obs ervations, which place 



estimates of D/H at 2.78 



+0.44 
-n as 



x 10" 



Kirkman et"ai]|2003ft 



10" 5 (lO'Meara et alj|2006h . However, de 



and 2.82+^ 

creasing the initial deuterium abundance in our models leads 
to less robust HD-cooling, and so only serves to strengthen 
our central conclusion that metal-free gas is unlikely to be 
cooled, by HD, to temperatures close to Tcmb- 

In §|4] we discuss recently updated rate coefficients for 
some of the most important reactions, and how their imple- 
mentation affects our results. 



2.2.3 Thermal Model 

The following processes are included in the net cooling rate: 
collisional excitation and ionization (of H, He, and He + ), re- 
combination (to H, He, and He + ), dielectric recombination 
(to He), Bremsstrahlung, Compton coolingjf] and molecular 
cooling by H2 and HD. In practice, the last two processes, 
as well as collisional excitation of HI, dominat e in our calcu- 
lation s. We adopt the expression provided bv lGalli fc Pallal 
(1998) for H2 cooling. In the fossil gas case, the HD-cooling 
rate is calculated using the analytic fit for low densities 
(n < 10 3 cm -3 ) given by equation (5) in iLipovka et al.l 
(2005). In the spherical collapse runs, we adopt the length- 
ier polynomial fit (equation 4 in the same source), which is 
accurate for gas densities n < 10 s cm -3 . 

We note that Tomb is a "temperature floor," below 



3 Num erical expressions for these cooling processes can be found 
in, e.g., lHaiman et~al1 lll996h . 



which gas cannot cool radiatively; if the gas temperature 
were below Tomb, interaction with photons in the roto- 
vibrational bands would heat, rather than cool the gas. In 
order to mimic this behavior, we multiply Ah 2 and Ahd by 
a correction factor (T — Tcmb) / (T + Tcmb)- This ensures 
that cooling is shut-off as the temperature approaches Tcmb 
from above (wh ereas the correction beco mes negligible when 
T S> Tcmb; see lHaiman et allll996l and I Johnson fc Bromml 
2006 for a somewhat more accurate approach). 

Our thermal model includes heating from photo- 
detachment of H~, high energy electrons resulting from 
photo-ionization of helium (see, e.g., Appendix B in 
lHaiman et al.l |l996), as wen as compressional heating in 
the model of adiabatic collapse (see OSH08 and references 
therein for more details). In practice, the latter dominates 
in the regime of our calculations. 

In order to follow the coupled chemical and thermal 
evolution of the gas we use the Livermore solver LSODAR 
to solve the stiff equations. 



2.3 HD and H2 Photo-dissociation in the 
Optically Thin Limit 

HD and H2 can be dissociated by photons with energies 
in the range 11.2-13.6 eV, to which the universe is largely 
transparent even before the IGM is reionized. Although both 
HD and H2 have LW lines above 13.6 eV, we do not include 
photons above this energy, because they will have been ab- 
sorbed by neutral hydrogen elsewhere in the IGM prior to 
reionization. Here we describe the details of the calculation 
for HD photo-dissociation; however, the calculation for H2 
is entirely analogous, so the following applies equally well to 
both molecules. 

Excitation of the HD molecule to its B 1 and C 1 H U 
electronic states and subsequent radiative decay leads to dis- 
sociation when the system decays to the vibrational contin- 
uum of the ground state, rather than back to a bound state. 
Here we discuss the optically thin case, in which the pro- 
cessing of the LW spectrum by HD itself (as well as by H2 
and HI) is assumed to be negligible. The dissociation rate 
for molecules initially in the electronic ground state with 
vibrational and rotational quantum numbers (v, J) is given 
by: 



' fdis 



(•>) 



where /dissV,./' is the dissociation probability from the ex- 
cited state («', J') and the pumping rate is given by 



J,v',J' 



' Jv 

4na v — — dv. 
hpv 



(6) 



Here a v is the frequency-dependent cross-section of a given 
transition, hp is Planck's constant, and the specific inten- 
sity just below 13.6 eV is hereafter normalized as J„ = 
J21 x [10 -21 erg cm -2 sr _1 s _1 Hz -1 ]. As mentioned above, 
in our model there is a sharp cut off in the radiation spec- 
trum above 13.6 eV; the lower limit, v t h, is the frequency 
threshold, corresponding to the longest-wavelength photons 
included, A~ 1105 A. 

In principle, the total dissociation rate also depends on 
the level populations of the molecule, which in turn depend 
on the incident radiation field as well as the temperature and 
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density of the gas; thus, the total dissociation rate should 
be: 



(7) 



where f Vj j is the fraction of molecules initially in the ro- 
vibrational state denoted by v, J. For simplicity, we assume 
that all HD and H2 molecules are initially in the ground 
state (i.e. f v=0tJ=0 = 1). This is a reasonable approximation 
for low gas densities, at which the populations of higher 
ro- vibrational states are very small. However, the level pop- 
ulations of H2 and HD reach their values at local thermody- 
namic equilibrium when gas densities rise to n > 10 4 cm -3 , 
and n > 10 6 cm -3 respectively. In §[3] we discuss the differ- 
ences in the dissociation rates if both molecules are assumed 
to be in LTE, and how this impacts the results discussed be- 
low. 

We include 28 discreet spectral lines of HD and 25 of 
H2 , all involving transitions from from the ground electronic 
state X 1 t° the B 1 and C 1 Tl^ excited states. We use 
the necessar y data for the Lyman an d Werner bands of HD 
provided by lAbgrall fe Rouefj J2006TI. For those of H 2 , the 
relevant data were t aken from Abgrall et al.l (|l993aT ). and 
lAbgrall et all (Il993bh. We use the updated dissociation frac- 
tions for H2 in lAbgrall et al.l |2000l ). The numerical wave- 
length resolution in the calculations (AA = 5.8 x 10 _5 A at 
the lowest temperatures) is sufficient to resolve the Voigt 
profile of each line and explicitly account for overlap of 
the Lorentz wings. We find the following photo-dissociation 
rates in the optically thin limit: fcdiss.HP = 1.55 x 10 9 J„ s _1 , 
and fc d i ss ,H 2 — 1-39 x 10 9 Jn s - 1 , in excellent ag r eemen t 
with those found previously by Idover fc Jappsenl l|2007f) : 
fcdi ss ,HD = 1.5 x 10 9 J s" 1 , and fc d i ss ,H 2 = 1.38 x 10 9 Ji; s" 1 . 
Here J„ denotes the intensity at the average LW band of 
HD and H2, with energy 12.4 eV, as discussed above. 



2.4 Self-Shielding of HD and H 2 

When sufficiently high column densities of HD or H2 
build up, (Nnr>,Nn 2 > 10 13 cm -2 ), the LW bands be- 
come optically thick and the rates of photo-dissociation 
are suppressed. We parameterize this effect by a shield 
factor, /shield, akin to that given by DB96 in their 
study of H2 self-shielding. In particular, /shield, hd = 

fcdis S ,HD(A r HD)Adiss,HD(A r HD = 0) where fc diss , HD (iVHD = 0) 

is the dissociation rate in the optically thin limit (equa- 
tion [7|, and the shield factor for H2, /shield, h 2 , is analo- 
gously defined. Our treatment of H2 self-shielding differs 
from the previous study by DB96 in that we assume all H2 
is in the ro- vibrational ground state (as described above), 
while the latter used a model allowing for populations in 
higher ro- vibrational levels due to collisional excitation and 
"UV pumping" by the incident radiation field. Nonetheless, 
we find that a good analytical fit for both H2 and HD self- 
shielding is provided by the same functional form as equa- 
tion (37) for /shield, h 2 m DB96. We also find that the self- 
shielding behavior of the two molecules is nearly identical, 
(see Figure [TJ, which might be expected on the basis of 



HD self — shielding 
U z self— shielding 
HD shielding fit 
H. shielding fit 



0.01 =- 




: o : -° 



Figure 1. Solid (blue) and broken (green) curves show the nu- 
merically calculated self-shielding factors for HD and H2 respec- 
tively; these are denned as the ratio of the dissociation rate at a 
given column density and the optically thin dissociation rate, and 
are shown as functions of the respective column densities. Dashed 
(orange) and dotted (magenta) curves show the values obtained 
by using the fitting formulae from equation [8] for HD and H2 
self-shielding, respectively (the slight difference here arises only 
because of the different masses of the two molecules, which mod- 
ifies the Doppler parameter). All are shown at T = 200 K. 



the similarity in their electronic structures Q; thus, we use 
the following fitting formula for both /shield, h 2 (A r H 2 , T) and 

/shield, HD (NhD,T): 



/shield (N,T) 



0.9379 



(l + x/Ds) 



0.03465 
(1 + x) - 47 



x exp [-2.293 x 10" 4 (1 + x) c 



(8) 



where x = AT/8.465 x 10 13 cm -2 , N is the column density 
of the self-shielding species, D5 = bo/10 5 cm s _1 , and the 
Doppler broadening parameter, bp, depends on the mass of 
the molecule (which accounts for the slight difference in the 
self-shielding formula for the two molecules), as well as the 
temperature. 



2.5 Shielding by HI and Mutual Shielding of H 2 
and HD 

In addition to self-shielding, HD and H2 can also shield each 
another, and both can be shielded by HI, which has absorp- 
tion lines in the range 11.2-13.6 eV; thus, suppression of 
the photo-dissociation rates depend on the relative strengths 
and positions of the HD, H2, and HI lines, as well as the col- 
umn densities of each species, A r HP, Nn 2 , and Nhv 



4 The similarities in the line strengths of H2 and HD, quantified 
by the product of the oscillator strength and dissociation fraction, 
can be seen in Figure [2] below. 
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Figure 2. The figure illustrates the wavelengths and strengths 
of the relevant LW lines of H2 and HD, as well as the Lyman 
series of HI. The solid (purple) and dashed (green) lines indicate 
the product of the oscillator strength and dissociation fraction for 
HD and H2 respectively, drawn at the position of each line center. 
The wavelengths and oscillator strengths of the hydrogen Lyman 
series are similarly shown by the blue dash-dot lines (a dissoci- 
ation fraction in this case is not applicable). At high HI column 
densities, the HD lines that dominate the dissociation rate are 
those at at ~ 960, ~ 980, and ~ 990 A, and the dominant contri- 
butions to the H2-dissociation rate are made by the absorption 
lines at ~ 945, ~ 980, and ~ 990 A. 



We include the first nine Lyman HI lines in the rel- 
evant wavelength range; while the line center of the Lya 
line is outside this range, we nevertheless include it, as its 
contribution to the shielding becomes important due to line 
broadening at high HI column densities. For illustration, in 
Figure [2] we show the positions and strengths of the most 
significant lines of each species. 

2.5.1 Shielding of HD by H 2 and HI 

Taking shielding into account, the contribution to the HD- 
dissociation rate for a particular line at frequency v becomes: 

fcdiss,HD,!/(iVHD) = fcdiss,HD, I /(A , HD = 0) eXp(— T„), (9) 

where the optical depth is given by 

t u = cthd.^hd + crm.vNm + an 2 ,uNn 2 - (10) 

While HD self-shielding becomes important for A^hd > 
10 13 cm -2 , the offsets in the wavelength of the neighboring 
absorption lines (typically of order ~A) prevent H2 and HI 
from effectively shielding HD until their column densities 
are very high, Nu 2 > 10 20 cm -2 , and N m > 10 23 cm -2 . At 
these critical densities, which are essentially independent of 
A?hd, the H2 lines start to significantly overlap and the opti- 
cal depth due to H2-shielding is ~ a few at all wavelengths. 
In Figure [3] we show the evolution of all three column den- 
sities, as a function of the particle number density, in our 



Figure 3. Column densities reached in the spherical collapse runs 
for halos with virial temperatures above and below 10 4 K, shown 
by solid and dotted lines respectively. Both runs assume no back- 
ground flux (J = 0). 



Table 1. Coefficients for the fitting formula for 
/shicld,HD(A r HD,A r H 2 .^ r Hl,T), representing the total com- 
bined shielding factor of HD, including shielding by H2 and HI 
(equations 1111 and 1126 . The analytic fit for self-shielding is given 
in equation [8] 



Species 


a 


P 


7 (cm 2 ) 


1. HI 


1.620 


0.149 


2.848 X 10 23 


2. H 2 


0.238 


0.00520 


2.339 x 10 19 



one-zone collapse runs with J = for halos with virial tem- 
peratures both above and below 10 4 K (see below). This fig- 
ure shows that all three column densities reach the values 
where shielding becomes efficient, and thus a full calculation 
of the combined self-shielding is warranted. 

We provide a fitting formula to model the total shielding 
of HD, /shield, hd = /shidd.HD^HD, Nn 2 , Nm,T) in Table 1 
and equations II II and [T2] below, which is accurate to within 
a factor of two over a wide range of column densities, i.e. up 
to Nhd ~ 10 20 cm -2 , A^h 2 ~ 10 22 cm -2 , Nm ~ 10 24 cm -2 
and gas temperatures up to w 10 3 K (HD-cooling is unim- 
portant at temperatures above this value in any case): 



/shicld.HD = /shield (JVhd, T) x /1 (Nm) x h {N H2 ) (11) 

1 



(l + xO c 



x exp(-^i xi) . 



(12) 



Here Xj = Ni/ji and the index i takes the value i = 1 or 
2 to denote the relevant quantity for HI or H2 respectively. 
The coefficients a, /3, and 7 are given in Table 1. 

This fit is described by a product of three sep- 
arate functions, / s hicid,HD(A r HD, T), /shieid,H 2 (A r H 2 ), and 
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/shield, hi(^Vhi), which represent the shielding of HD due to 
each of the three species alone (eq. QT] above). Note that 
since the H2 and HI lines shield HD by their Lorentz wings, 
rather than their thermal cores, these factors (unlike HD 
self-shielding) do not depend on temperature. 

In general, one does not expect that the combined 
shielding factor is separable into a simple product of the 
three individual shielding factors. The full expression for 
/shieid,HD(-/VHD, ATh2, Nm, T) is a sum over all of the individ- 
ual cross-sections of the HD lines, each suppressed by the 
frequency-dependent total optical depth (eq. I10[) . divided 
by the optically thin rate. However, we have found that in 
practice, when suppression by nearby HI lines is negligible 
(iVm < 10 23 cm -2 ), one HD line is much stronger than 
all others (at w 950A, see Figure [5]). Since this single line 
dominates the dissociation rate, the shield factor reduces to 
the simple product (eq. [11] above). In the regime of rela- 
tively strong HI shielding (Nm ~ 10 24 cm -2 ), a few HD 
Lyman lines together dominate the dissociation rate. How- 
ever, we find that the product a\ x /diss, a for these lines (at 
~960, ~980, and ~990 A) are similar. If we approximate 
that these lines have identical strengths, the total shielding 
factor again reduces to the simple product in equation 1111 
above. Because these line strengths are not precisely equal, 
the largest discrepancies in the product formula and 'true' 
shielding behavior are seen at TVhi = 10 24 cm -2 . However, 
in general, we find that this simple product is accurate, to 
within a factor of ~two, at the low temperatures (T < 200 
K) and the high column densities of interest. 

Figure [4] shows the results of the exact shielding factor 
calculations for a gas temperature of T — 200K, and com- 
pares these to the analytical fits for a number of combina- 
tions of the three column densities. The largest deviations 
are seen in the bottom panel, for Nm = 10 24 cm -2 , and 
at H2 and HD column densities of JVh 2 = 10 22 cm -2 and 
10 14 cm -2 < Nhd < 10 15 ' 5 cm -2 . In this regime, the accu- 
racy of the fitting formulae is somewhat worse than a factor 
of two. However, in practice, this column density combina- 
tion - with relatively low Nho and exceedingly high values 
of both Nh 2 and Nm - does not occur in our calculations 
(see Figure [3]). 



2.5.2 Shielding of H 2 by HI and HD 

The shielding of H2 by HD and HI is entirely analogous 
to that discussed in the preceding section, so we will limit 
the discussion here to a few noteworthy points. Most im- 
portantly, we find that HI shielding of H2 is nearly identical 
to HI shielding of HD; accordingly, we model both with the 
fitting formula for / 8 hield,Hl(^VHi), given by equation [12] and 
Table 1. The explanation for this is similar to that given 
above; namely, the relative positions of the H2 and HI lines 
is such that only the wings of the HI lines shield H2 when the 
column densities of Nm are sufficiently large (Nm > 10 23 ). 
Because the H2 and HD lines are comparably spaced rela- 
tive to the HI lines, the shielding effect of HI should indeed 
be similar for both. 

When the HI column is below the critical level for 
strong shielding of H2, we find that a few Lyman lines to- 
gether dominate the dissociation rate, and that the product 
o~x x /diss, a for these lines (at ~945, ~980, and ~990 A) are 
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Figure 4. The combined HD shielding factor, including self- 
shielding and shielding by H2 and HI. Several combinations of 
column densities are shown, as labeled, near the critical column 
densities for HI and H2 shielding. The solid curves show the ex- 
act numerical calculations, and the dotted curves show the values 
obtained from a fitting formula (equations 151 fTTT and [12] and Ta- 
ble 1). 
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Figure 5. The combined H2 shielding factor, including self- 
shielding and shielding by HI. Several combinations of column 
densities are shown, as labeled, near the critical column densities 
for HI shielding. The solid curves show the exact numerical cal- 
culations, and the dotted curves show the values obtained from a 
fitting formula (equations 181121 1131 and Table 1). 

similar. At larger neutral hydrogen columns (Nm ~ 10 25 ), 
a single Lyman H2 line makes the dominant contribution - 
by a large margin over all others - to the dissociation rate. 
Thus, the total shielding factor can again be simply modeled 
by a product of the shielding formulae, given in equations [8j 
[T3I and Table 1: 

/ sh icld,H2 = /shield (Nv 2 , T) X /1 (JVhi) (13) 

The results of the exact shielding factor and comparison to 
this analytical fit are shown in Figure[5]for a gas temperature 
of T = 200K. 

In principle, HD can also shield H2, but in practice this 
effect will likely always be negligible, as the HD column den- 
sity is typically dwarfed by those of H2 and HI0. Thus, our 
treatment does not include HD shielding of H2. 



3 RESULTS 

To assess whether HD-cooling can be suppressed by a persis- 
tent LW background, we ran our one-zone models at various 
different specific intensities J21 . Unless stated otherwise, the 
spectrum of the radiation is modeled as a black-body with 
a temperature of 10 5 K, approximating the hard spectrum 
expe c ted to characterize Pop HIT stars ( Tumlinson &: Shull 
120001 : iBromm et aT1l200ll ; ISchaered I2OO2V For comparison, 
in § 13.2.21 below, we investigate the effects of illumination 
by a cooler blackbody, T ~ 10 4 K, intended to represent 
the softer spectrum of a more typical metal-enriched stellar 
population. These are referred to hereafter as types 'T5' and 
'T4' respectively (SBH10). 

5 The fractional abundance of HD relative to that of H2 could ex- 
ceed the cosmological D/H ratio by a large factor, o wing to chem- 
ical f ractionation at low temperatures (see, e.g., iGalli fc Pallal 
1998). However, it never exceeds ~ 10 in our models. 



We use a Newton-Raphson scheme to determine the 
strength of the LW radiation required to keep the gas 
temperature greater than a factor of ~ 2 above that 
reached in the absence of any LW radiation, T m i n ,j=o, on 
the timescales described below; this is referred to here- 
after as the critical intensity: J cr it,HD (in the usual units 
of HT 21 erg cm" 2 sr" 1 s" 1 Hz" 1 ). 

3.1 HD-Cooling in Constant-Density Fossil HII 
Gas 

The first of the physical scenarios we consider is a "fossil" 
HII region, which could occur in a patch of the low-density 
IGM that has been photo -ionized and heated by a short- 
lived massive star (see, e.g. lOh fc H aiman 2003), or possibly 
in a denser shell of primordial gas, compressed by shocks 
from a supernova (SN). 

We are interested in whether gas with such fossil ioniza- 
tion can cool efficiently in the presence of a LW background 
and return, via HD cooling, to a state close to its initial low- 
entropy state, prior to the ignition of the ionizing source (or 
prior to its shock heating). 

In a first set of runs, we assume that the number den- 
sity remains constant at the low value of n = 10~ 2 cm" 3 , 
characteristic of a slightly over-dense (by ~ a factor of 10) 
region of the IGM at z = 10. We find that in such a rarefied 
patch, the gas is not able to cool on a realistic timescale, 
because the HD-cooling time is longer than the present age 
of the universe even in the absence of any LW background. 
The thermal evolution of the gas in this case is shown by the 
right set of (purple) curves in the upper panel of Figure [5] 
The corresponding fractional abundances of electrons, H2, 
and HD are shown in the lower panel of the same figure. 
The figure extends to a total elapsed time that exceeds the 
Hubble time, and shows that there is, technically, a critical 
flux (Jcrit.HD ~ 10" 6 ) that would suppress the HD-cooling 
that would otherwise occur after a few x 10 18 seconds. This, 
of course, is unphysical, and in practice, the question of 
HD-cooling being suppressed by UV radiation is moot for 
such low-density gas. Nevertheless, the figure illustrates the 
chemical/thermal behavior, and also provides a useful check 
on our code (see below). 

In the second set of runs, the total particle number den- 
sity was set to a higher value of n — 10 2 cm -3 . This is an 
unphysically high density for a characteristic "flash-ionized" 
fossil region in the low-density IGM, but may represent pri- 
mordial gas compressed by SN shocks. In the no-flux case , 
we reproduce the main result of IJohnson fc Bromml (2006), 
namely that HD-cooling allows the gas to reach the temper- 
ature of the CMB in a time that is shorter than the Hubble 
time. This gas, however, is optically thin to radiation in the 
LW bands, and dissociation of HD (as well as of H2) is ef- 
ficient. We find that for J21 > 10~ 2 , the gas cannot cool to 
temperatures less than ~ 200 K. This is shown in the top 
panel of Figure [6] by the left set of (blue) curves. 

The above two cases (n — 10 -2 and n — 10 2 cm -3 ) 
serve to illustrate an important point (and a check on our 
code). All of the relevant timescales for the system, includ- 
ing the formation timescale for HD (tform), the HD cool- 
ing time (f coo i), and the (HH+e — s-HI) recombination time 
(f rC c), scale as 1/n. The exception is the photo-dissociation 
timescale, which scales with the flux strength, idiss oc JJ[ , 



© 2010 RAS, MNRAS 000, 000-000 



8 J. 



Wolcott- Green and Z. Haiman 



2 2 



Q. 5 



-10 



_ 1 1 . 


— |— i i | i 




wi H -10* 


\n H -10 • 









1 3 


- T 

1 CMB 

J 21 -0 

10 3 ' 


j.o- 2 






\. 10"' - 
J al =0~ ^ 


- — 1 — — 1 — 1 — \ 




H 1 1 1 ^ 

n„-10- a 


e 








0\ \ : 


H 2 


/ . 




HD 

- ■ 1 , , 


,1,1, 


i i i i i _ 


12 


14 16 
log(t/ sec) 


18 20 



Figure 6. The top panel shows the thermal evolution of constant 
density gas that is initially fully ionized. The right/left set of solid 
purple/blue curves adopt densities of 10 — 2 and 10 2 cm -3 , appro- 
priate for a flash— ionized patch of the IGM and for primordial 
gas compressed by a SN-shock, respectively. In both cases, we 
show the evolution when the gas is exposed to LW backgrounds 
of various intensities, as labeled. The dashed lines show the tem- 
perature reached in deuterium- free gas, in the absence of any LW. 
The temperature of the CMB, Tomb = 2.7(1 + z), is shown by 
the dotted blue curve (starting from z = 20). The bottom panel 
shows the fractional abundances of electrons (light blue dotted 
lines), HD, and H2 (solid blue and green lines respectively) in 
the low— density case. The fractional abundances are not shown 
in the case of the shock-compressed gas for the sake of clarity, 
and because the patterns they follow are identical to those shown 
in the lower panel. This figures shows that there is a critical flux, 
^crit,HD = 3.8 X 10 — 2 (n/10 2 cm -3 ), that suppresses HD-cooling 
and prevents the gas from reaching the CMB temperature. 

in the absence of shielding. Thus, the history of the system 
should only depend on J2i/n when the time is rescaled ac- 
cordingly^ This simple scaling is evident by the two sets 
of (purple and blue) curves in the top panel of Figure [6] 
Most importantly, there is indeed a critical flux that pre- 
vents the gas from reaching the CMB temperature by HD 
cooling; in constant density gas with a high initial elec- 
tron fraction, we find the value of this flux is J cr i t ,HD = 
3.8 x 10~ 3 (n/10 2 cm" 3 ). 

Finally, an interesting question is whether HD-cooling 
prevented, for the cases in which J21 exceeds the critical 
value, by direct photo-dissociation of HD, or the inabil- 
ity of sufficient HD to form due to H2 photo-dissociation. 
To answer this question, we performed runs in which the 
H2-dissociation was artificially turned off. In these runs, 
we find that the gas is still able to cool to Tcmb for 

6 The interested reader can find a much more det ailed discussion 
of this point for the analogous case of H2 cooling in lOh fc Haimar] 
120021 . 



J21 ~ 4 x 1CP 2 , illustrating that the LW flux prevents 
HD-cooling via Hi destruction, rather than via direct HD- 
dissociatio n. This point has been discus sed by previous au- 
thors, e.g. lNakamura fc Umemural |2002f ) showed that a crit- 
ical abundance of H2, xn 2 > 10~ 3 is required for the gas to 
reach sufficiently low temperatures for HD to become the 
dominant coolant (T < 150 K), and therefore H 2 dissocia- 
tion can prevent HD-cooling l|Yoshida et al.| [2007'l . The min- 
imum temperature reached by fossil ionized primordial gas 
in the absence of HD, and the dependence of this minimum 
temperatur e on gas density and L W flux was also discussed 
in detail bv lOh fc Haiman] (|2003T ). 



3.2 HD Cooling in Collapsing Halos 

3.2.1 HD Cooling in Halos with T vir > 1Q A K 

It has been shown that primordial gas in the late stages 
of runaway gravitational collapse can reach temperatures 
close to Tcmb via HD-cooling , provided that a large in itial 
ioni zation fraction exists (e.g. IJohnson fc Brommll2006l ; see 
also iMachida et~aH I2005T I . 

This scenario can be realized in sufficiently massive 
halos, which are collisionally ionized upon shock-heating 
to their virial temperatures. The post-shock gas can cool 
faster than it recombines, leaving a large out-of-equilibrium 
elect ron fraction to catalyze both H2 and HD formation 
(e.g.lShapiro fc Kandll987l ; ISusa et al.lll998l ; lOh fc Haimanl 
l2002h . It may also be the case that "pre-ionized" halos ex- 
ist within fossil HII regions, which wi ll undergo a phase 
of ef ficient HD-cooling upon collapse l|johnson fc Bromml 
2006). This scenario, however, is less plausible: a halo large 
enough to remain bound once photo-heated (to T > 10 4 K) 
may be difficult to completely ionize, as the large HI col- 
umn densities will le ad to non-negligible HI self-shielding 
|Diikstra et all 12004 ); flash-ionization by a single short- 
lived star (required to allow subsequent recombination and 
cooling) is even less likely. 

Regardless of the nature of the initial ionization, Fig- 
ure [7] shows the thermal evolution of the collapsing gas ex- 
posed to LW backgrounds of various intensities. The initial 
number density is set to the characteristic baryon density in 
halos upon virialization, 

n^0.3cm- 3 (i±^) 3 , (14) 

and the gas begins cooling from the temperature T ~ 10 4 K 
(quickly established either by a period of photo-heating, or 
by shock-heating to near the virial temperature, accompa- 
nied by rapid HI cooling). 

Of primary importance in this case - as opposed to 
the fossil gas discussed in the previous section - are the 
large column densities of HD and H2 that build up and 
shield both populations against the LW background. Fur- 
thermore, the collapse itself leads to more efficient forma- 
tion of both molecules because the formation timescale, as 
mentioned above, scales as tform °c 1/n. Consequently, the 
critical intensity should be larger than that found for the 
low-density fossil gas. This is indeed borne out by our re- 
sults; nevertheless, as the comparison between the solid and 
the dashed curves (the latter representing deuterium-free 
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Figure 7. Thermal evolution of initially ionized gas, collapsing 
in a massive halo (T v j r > 10 4 K) exposed to LW backgrounds 
in the range near the critical value, J cr j t hd = 3.6 X 10 — 1 (solid 
curves). The incident spectrum is that of a blackbody with a tem- 
perature of 10 5 K, characteristic of massive Pop III. 1 stars. The 
turnaround rcdshift is set to z = 17 and the temperature is ini- 
tialized T 10 4 K. The dashed curves show, for comparison, the 
temperature evolution in deutcrium-free gas, exposed to the same 
LW fluxes. The dotted curves show the thermal evolution when 
HD-dissociation is artificially switched off, for the same values of 
J21, as labeled. This illustrates that, except in a small range of 
flux intensities, the destruction of H2 by LW photons, rather than 
direct HD photo-dissociation, is the primary factor determining 
the minimum temperature reached in the gas. The temperature 
of the CMB is shown by the blue dotted curve. 



gas) in Figure [7] shows, the effect of HD-cooling is still al- 
most entirely erased for the relatively low values of J > 1. 
The critical intensity in this CcLSG, ELS defined above, is found 
to be Jcrit = 3.6 x KT 1 . 

This threshold value is most notable for being approxi- 
mately five orders of magnitude lower than the critical flux 
required to completely suppress H^-cooling in the same ha- 
los. As shown in SBH10, the latter critical flux in halos with 
T v i r > 10 4 K is Jcrit, h 2 > 10 . This large critical flux corre- 
sponds to the value that results in an H2-photo-dissociation 
rate that matches the H2 formation rate, at the critical den- 
sity of n ~ 10 4 cm -3 of H2 (see the earlier work by O01 
for a detailed discussion of the physics determining J cr it,H 2 
in primordial gas without ionization/shock-heating). This 
Jcrit, h 2 ~ 10 4 separates gas in which H2-cooling is fully sup- 
pressed (with the gas temperature remaining near ~ 10 4 K) 
and halos in which H2 cooling significantly lowers the tem- 
perature. A point that was also found (but not emphasized) 
by SBH10 (and also by O01) is that even for J21 well be- 
low Jcrit, h 2 , the minimum temperature to which the gas can 
cool via H2 can be significantly elevated. This is also clearly 
visible in the deuterium-free runs in Figure [7] the mini- 



mum temperature is ~ 150K for J21 = 0, but is elevated to 
~ 300K already for J21 = 1. 

The behaviour of the gas, and the reason for the 
elevated temperature, can be described as follows (see 
a detailed discussi on in the constant-density case in 
lOh fc Haimanll2002f ). Starting from T ~ 10 4 K, the gas ini- 
tially cools via H2 and recombines on time-scales much 
shorter than either the photo-dissociation or the free-fall 
timescale. However, when the temperature is lowered to a 
J-dependent critical value of a fewxl0 3 K, H2-dissociation 
becomes important, limiting the H2 abundance, and reduc- 
ing the cooling. The cooling time eventually becomes compa- 
rable to the free-fall time, resulting in the sharp turn away 
from the nearly vertical directions of the n — T curves at 
the initial density in Figure [7] For higher fluxes, this subse- 
quently results in an elevated gas temperature (at fixed den- 
sity). Eventually, the compressional heating rate becomes 
equal to the H2-cooling rate, setting the temperature min- 
imum. It is worth noting that for sub-critical values of 
J21 < Jcrit, hd, HD is able to cool the gas to temperatures 
below 150K, but the LW background still has the subtle ef- 
fect of raising the minimum temperature to which the gas 
can cool via HD. 

As in the constant density case, an interesting question 
is whether ultimately the HD-cooling is controlled by direct 
photo-dissociation of HD, or by H2-dissociation. To answer 
this question, we repeated the runs shown in Figure [7] under 
the same conditions, but with HD dissociation artificially 
switched off (the results are shown by the dotted lines in 
Figure [7)). In this case, we find again that (except in a nar- 
row range of LW intensities near J21 ~ 4 x 10 _1 ), direct 
photo-dissociation of HD does not determine the minimum 
temperature reached in the gas. Rather, it is the diminished 
abundance of H2 in the presence of the LW background that 
regulates the abundance and thereby the cooling efficiency 
of HD. 

We have also investigated the thermal evolution of the 
gas when H2 is artificially prevented from dissociating (not 
shown), but the physical set-up is otherwise analogous to 
the runs (in which deuterium is included) shown in Figure [7) 
This "academic" exercise is useful in order to determine the 
critical flux that would prevent HD-cooling entirely by di- 
rect HD photo-dissociation. In the analogous case for H2, as 
mentioned above, J cr it,H 2 is traditionally defined as the spe- 
cific intensity capable of completely suppressing H2-cooling, 
thereby preventing the gas from falling below the tempera- 
tures reached by atomic line cooling, T ~ 8 x 10 3 K (O01, 
OSH08, SBH10). We find that for J > 6 x 10 4 , the cooling 
history of the gas is nearly identical to that of deuterium- 
free gas in the absence of a LW background. Thus, the 
intensity required to fully suppress HD-cooling by direct 
HD-dissociation is comparable to that of H2 (the latter was 
found by the latest studies (SBH10) to be ~ 1.2 x 10 4 , but 
we find a factor of ~ five greater critical value for H2; see 
below). In fact, this is not surprising, given that the photo- 
dissociation timescales (tdiss = [kdiss{N = 0) x /shield] -1 ) are 
very similar for the two molecules. 



3.2.2 Tvir > 10 4 K Halos Illuminated by 'T4 ' Radiation 

Up to this point, we have considered only incident radiation 
with the hard spectrum expected to characterize Pop III. 1 
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Figure 8. Thermal evolution of gas collapsing in a massive halo, 
exposed to different LW backgrounds, as in Figure [7] except the 
incident spectrum is that of a softer blackbody with tempera- 
ture T = 10 4 K, representing a more typical metal-enriched stel- 
lar population. This softer spectrum contains many more photons 
down to the photo-dissociation threshold of H — at 0.76eV. The 
enhanced rate of H~ photo-dissociation reduces the critical flux 
that prevents HD-cooling by more than an order of magnitude 
compared to the harder spectrum, to J C rit,HD = 10~ 2 . 

stars ('T5'). In the case of a collapsing halo, it is reasonable 
to ask about the effects of irradiation by the more the typical 
stellar spectrum ('T4'), on HD-cooling. 

Figure [8] shows the temperature evolution in the gas ir- 
radiated by a T4 spectrum of various intensities. It is clear 
that HD is considerably more fragile in the presence of this 
softer spectrum, withstanding a LW flux no greater than a 
feeble J cr it,HD = 1CP 2 . This is not surprising in the light of 
studies that have found the same effect for H2 (O01; OSH08; 
SBH10), namely, that H2-cooling is much more effectively 
suppressed by the T4 type spectrum. The reason is the di- 
minished abundance of hydride (H~), an intermediary in 
the formation of both H2 and HD. Hydride, whose ioniza- 
tion threshold is 0.76eV, is more efficiently photo-dissociated 
by the softer spectral type (O01, SBH10). This again is a 
manifestation of how an external radiation field can regu- 
late HD-cooling via the destruction of an intermediate in its 
formation pathway. 

3.2.3 HD Cooling in T vir < W 4 K Halos 

It has long been known that pristine gas in the first mini- 
halos canno t form sufficient H2 to cool below a few hundred 
Kelvin (e.g. lHaiman et aT1ll996M . Because free electrons act 
to catalyze the formation of HD, as well as H2, it is not sur- 
prising that HD abundances remain to o low in such halos 
to pl ay a significant role in cooling (e.g. Ijohnson fc Bromml 
120061 . and references therein). This is borne out by our re- 



sults, shown in Figure [9l for a halo that begins collapsing 
from a temperature of T « 20K at the turnaround redshift, 
2 = 17, and is illuminated by a 'T5' spectrum. This is the 
same configuration as in Figure [7] except that the initial 
temperature is assumed to be low (i.e., lacking any strong 
shocks able to collisionally excite or ionize the gas). 

We include a discussion of this scenario for the purpose 
of highlighting a few noteworthy points. First, we find a 
critical flux for full H2-dissociation in this case (as it is tra- 
ditionally defined, see § 13.2.11 above) of J cr it,H 2 = 6.1 x 10 4 , 
which is a factor of ~ five greater than that found in the re- 
cent study by SBH10. This difference owes to our use of 
the analytic fit for H2 self-shielding (equation [S|, which 
assumes all molecules are initially the ground state, while 
SBH10 used the shield factor fit (equation 36) from DB96, 
which gives a very good approximation to the self-shielding 
when the H2 roto-vibrational levels reach LTE populations. 
This illustrates an important point that will be discussed in 
greater depth in ijl] namely, self-shielding is less effective 
when higher ro-vibrational states of the molecule are pop- 
ulated (see more discussion in §[4] below). It is also worth 
noting that we find J cr it,H 2 — 6 x 10 3 - a factor of 2 lower 
than that found by SBH10 - if the self-shielding is mod- 
eled instead with the more accurate formula provided by 
DB96 (equation 37), rather than the power-law fit (equa- 
tion 36) used by SBH10. In general, the "real" LW intensity 
required to kill H2-cooling entirely will depend on the de- 
tailed ro-vibrational distribution of H2 molecules. The values 
of J cr it,H 2 found here using the ground-state vs LTE shield- 
ing treatments serve to bracket the range for the true critical 
threshold, with the former approximating the upper limit on 
Jcrit,H 2 - The threshold flux in this case is of particular inter- 
est as it has important implications for the number of halos 
that remain H2-poor, and thus for the abundance of halos 
collapsing directly to massive seed black holes, because such 
halos probe the rare, exponential tail of the fluctuating UV 
background (|Diikstra et al.|[200St ). (The interested reader is 
encouraged to see SBH10 for a detailed discussion of these 
issues.) 

Second, while there is a clear bifurcation in the cooling 
history of the gas around the J cr it,H 2 discussed above, we 
note that it does not cool to the low temperatures required 
for HD-cooling to become important (T ~ 150K) even when 
the intensity is well below this threshold (indeed, even in the 
absence of any LW background). Thus, again, H2 abundance 
plays the primary role in regulating HD-cooling, though in 
this case the result is simpler: HD-cooling never becomes 
important because H2-cooling is never strong enough for 
the gas temperature to fall below ~ 200K. 

Finally, in this case HI shielding serves to increase the 
threshold LW flux by ~ 40% above that which H2 could 
otherwise withstand. By contrast, J cr it,HD found in the two 
preceding sections is not sensitive to the effect of H2 shield- 
ing by HI, in spite of its strong dependence on the H2 abun- 
dance. This is because, in the models of T vir > 10 4 K halos, 
the HI column density does not reach the critical value above 
which it effectively shields H2 (iVm > 10 23 cm -2 ) until very 
late in the collapse (i.e., once the particle density has reached 
n > 10 5 ' 5 cm -3 ). However, sufficiently high column densities 
of HI do build up earlier in the high- J21 runs of TVi r < 10 4 K 
halos, resulting in the modest increase in J cr it,H 2 as quoted 
above. 
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Figure 9. Thermal evolution of gas collapsing in a halo exposed 
to LW radiation (with a T5 spectrum) in Figure [7] except that 
the initial gas temperature is assumed to be much lower, 2 IK. 
This is relevant to minihalos (T v ; r < 10 4 K), or to the case when 
there are no strong shocks able to collisionally excite or ionize the 
gas. 

The results for the halo illuminated by a T4 spectrum 
are omitted here because they do not significantly differ 
from those described above for the T5 case, except that 
HVcooling is disabled by a much lower LW flux, as shown 
already by O01. 



4 DISCUSSION 

The most significant result found in this paper is that a UV 
flux can prevent HD-cooling from lowering the gas temper- 
ature to near that of the CMB. In particular the threshold 
value in collapsing halos, J cr it,HD ~ 4 x 1CP 1 , is approxi- 
mately five orders of magnitude lower than the critical flux 
required to completely suppress H2-cooling in the same ha- 
los. As explained in § 13.2.11 this large difference arises be- 
cause even for J21 well below J cr it,H 2 , the minimum tem- 
perature to which the gas can cool via H2 is significantly 
elevated, so that HD formation and cooling is not activated. 
What are the cosmological implications if a LW background 
prevents halos from cooling via HD line emission? The crit- 
ical flux we find should be compared to the level of the LW 
background Jb g expected at the redshifts we consider. For 
reference, an estimate of the background as a function of 
redshift is provided by requiring the number of UV photons 
produced by stars to b e sufficient to reionize the IGM. Such 
an estimate gives (e.g. iBromm fe Loebl 2003) 

Here the number of UV photons required to ionize hydro- 
gen is commonly taken to be iV 7 ~ 10, and the escape frac- 



tion / esc is the fraction of ionizing photons (13.6eV) escap- 
ing from galaxi es at high redshifts, whic h is expected to be 
close to unity (W halen et alj |2004| ; see iFernandez fc ShuH 
2010 for a recent discussion of the relevant / csc and for ref- 
erences to earlier works). From this equation, we find that 
the mean UV background at the time of reionization exceeds 
the critical flux J cl -it,HD by nearly an order of magnitude. 
While the radiation background will inevitably have spatial 
fluctuations, this implies that most halos collapsing in the 
early IGM, prior to reionization at z ~ 10 — 20, would be ex- 
posed to a super-critical flux and thus not able to cool below 
T ~ 200K. As a result, the emergence of low-mass PopIII.2 
stars (or stars comparable in mass to those formed in the 
low-redshift universe) would be postponed until supernovae 
polluted the IGM with heavy elements, and metal-line cool- 
ing subsequently enabled gas clouds to reach temperatures 
near Tomb- 

There are several issues that could be important for the 
above conclusion, which we have glossed over in the dis- 
cussion of our models for molecular shielding and one-zone 
spherical collapse. We next discuss some of these. 
Uncertainty in the Column Density of the Collapsing Re- 
gion. In our calculations, we have taken the diameter of the 
collapsing region to be of the order of the Jeans length. 
However, the effective size and column density will depend 
sensitively on the dynamical properties of the system, in- 
cluding bulk motions and internal velocity gradients in the 
gas, departures from spherical symmetry, etc. 

In order to address this uncertainty quantitatively, we 
have performed two additional sets of runs for a halo with 
T v i r > 10 4 K (analogous to those in § I3.2.1[1 . In the first, we 
increase the assumed size of the collapsing region by a factor 
of ten (i.e. to ten times the Jeans length Aj). This increases 
the column densities by the same factor, and accordingly, 
J cr it is larger by a factor of ~ 3 than the original result, 
due to more efficient self-shielding of H2 and - to a lesser 
extent - HD. Next, we investigated the cooling properties 
of a smaller collapsing core, with the assumed size reduced 
by a factor of ten, to O.lAj. In this case, a new effect arises. 
Namely, for the values of J21 at which the gas can (just) cool 
to around T ~ 200K, we find that direct HD dissociation is 
the dominant factor regulating the minimum temperature 
ultimately reached. In particular, switching HD dissociation 
off by hand in these cases decreases the gas temperature by 
a factor of ~ 1.5 at high number densities (n > 10 4 cm -3 ) 
for J = 10 _1 . When the flux is weaker than this, artificially 
disabling the HD-dissociation has little effect; for these low 
fluxes, however, the suppression of HD-cooling is modest to 
begin with. As expected, the critical flux is decreased in this 
case because the smaller column densities leave both H2 and 
HD more susceptible to dissociation. In particular, we find 
a critical value of J cr it ~ 10 _1 . 

The Impact of 3-D Gas Dynamics on Self-Shielding. It is 
important to note that a full treatment of the three- dimen- 
sional dynamics of the system and the complexities inherent 
in radiative transfer is needed to solve the shielding problem 
exactly. Our calculation is based on a model of a uniform 
slab of gas with no internal velocity (or temperature) gradi- 
ents. This is likely to be a poor approximation for a region 
undergoing runaway gravitational collapse, in which high gas 
velocities can produce significant Doppler shifts of the LW 
absorption bands of H2 and HD. In general, we expect this to 
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reduce the effective column densities, and the importance of 
shielding, compared to our calculations. This correction can 
be mitigated by the broadening of absorption lines at high 
column densities, which leads to line widths that are much 
larger than the average Doppler shift. In general, however, 
taking into account the possibility of Doppler shifts leads 
to the conclusion that our self-shielding results, and in turn 
the values quoted for J cr it,HD, are upper limits for both. This 
strengthens our argument above, namely, that most collaps- 
ing halos will see a super-critical flux. 

Resonant Scattering of Incident LW Photons. An additional 
effect ignored by our self-shielding calculation is that after 
absorbing a LW photon, a fraction of molecules will decay 
directly back to the original ground state, as opposed to cas- 
cading through a series of lower energy decays as we implic- 
itly assume. These photons are thus not eliminated, making 
the background flux s tronger than we calculate. However, 
Idover fc Brandl l|200ll ) estimate that in the case of Eh, such 
resonant scattering constitutes only a small fraction, 4-8%, 
of all LW absorption events (depending on the initial level 
populations; in particular, on the ortho/para ratio). Hence 
this is a minor effect, which again makes our conclusions 
about the suppression of HD-cooling conservative. 
Molecular Level Populations: Implications for Photo- 
dissociation Rates and the Critical Flux. As described in 
H2.3I our fiducial self-shielding calculations for HD and H2 
assume that all molecules are initially in the ro-vibrational 
and electronic ground states. In reality, molecules will oc- 
cupy higher ro-vibrational states due to collisional excita- 
tion, and this can significantly increase the rates of photo- 
dissociation. In particular, wc find that self shielding is less 
effective (i.e. the shield factor is larger by a factor of a ~ 
few) if we repeat our calculations, assuming LTE popula- 
tion distributions in the rotational levels (J 7^ 0) within the 
ground electronic and vibrational s tate. (Data for the HD 
and H2 energy le vels were taken from lAbgrall et al ]| 19821 and 
iDabro wski 1984] respectively.) Other studies have also noted 
that populations in higher vibrational levels (v 7^ 0) can 
significantly increase th e rates of photo-dissociation (e.g. 
Idover fc Jappsenl 120071 . and references therein). Thus, if 
the effects of collisional excitation are taken into account, 
Jcrit, hd/h 2 could be significantly lower than found above. 
The ro-vibrational distribution will be somewhat better ap- 
proximated by the ground state model up to the critical 
densities (n > 10 4 cm~ 3 for H 2 and n > 10 6 cnT 3 for HD). 
In our case, the fate of the collapsing clouds - whether it will 
ultimately cool to temperatures low enough for HD-cooling 
to become significant - is determined in the regime of some- 
what lower gas densities, below those at which equilibrium 
H2 populations are established, so that J cr i t values are likely 
closer to the upper end of the range. 

Uncertainties in the Gas Phase Chemistry. The preceding 
discussion has focused on various uncertainties associated 
with photo-dissociation rates of H2 and HD; however, the 
accuracy of any estimates of J cr it,HD can only be as good 
as the accuracy of the underlying chemical rate coefficients. 
Considerable attention has been dedicated to uncertainties 
in both H2 and HD chemistry and cooling, (e.g. |Savin et al.l 
l2004lGlover et al.ll200d lGloverll200sl , iGlover fe Abei .2008): 



Two reaction rates that are crucial for determining the 
abundance of H2, particularly in gas with a large initial ion- 
ization fraction, have been esti mated to be uncertai n by up 
to an order of mag nitude fe.g. lGlover fc Abejlioosh . These 
are the associative detachment channel for H2 formation, 
and the mutual neutralization of hydride and protons (reac- 
tions 10 and 13 respectively in the compilation by SBH10): 



H + H~ 



H + +H~ 



H 2 +e~ 



H + H. 



(16) 
(17) 



Both have been revisited recently and new the rmal rate coef- 
ficient s have been provided for rea ction (10) bv | Kreckel et al.l 
l^Old ) and for reaction (13) by IStenrup et all (|2009l V As 
noted above, the formation of HD occurs primarily via the 
reaction pathway shown in equations \T\- [3] so in general, 
the fractional abundance of HD is proportional to that of 
H2, justifying the emphasis here on uncertainties associated 
with the formation of H2. In order to quantify how these re- 
cently updated rate coefficients impact our results, we have 
performed additional runs for each of the physical scenar- 
ios discussed in § [3] The uncertainties these introduce into 
estimates of Jcrit, hd and the minimum temperature, T m i n , 
reached when J21 = 0, are summarized below. 

We examine three published rate coefficients for the 
mutual neutralization reaction; the largest of these at all 
temperatures is that provided by GP98, hereafter rCi 3 b. The 
rate used in our fiducial model (^13), originally provided by 
iDalgarno fc Leppl (|l987), is smaller than the others by up 
to ~ an order of magnitude for T > 10 3 K, and by a factor 
of several at lower temperatures . Q The value of the n ewest 
thermal rate coefficient, given bv lStenrup et all l|2009h . fci 3c , 
lies between those of fci 3 and fci 3b at all temperatures in the 
regime we study. Using fci3b, we find the critical flux, Jcrit, hd 
is a factor of ~ 2 lower than was found in each of the physical 
scenarios we have studied (see § [3] note that this does not 
apply to the model of halos with T v i r < 10 4 K), and T m i n is 
~ 30% greater in the spherical collapse models. This is easily 
explained: when a large ionization fraction exists, reaction 
(13) competes with reaction (10) for the common reactant, 
H~; thus, adopting a larger rate coefficient for reaction (13) 
leads to diminished abundance of, and less robust cooling by 
H2. Using the newest rate, fei3 C , we find the critical flux is 
~ 30 — 40% smaller than its original value for each physical 
scenario, and the minimum temperature is elevated ~ 20% 
above that found previously in the spherical collapse models. 
(Note that the minimum temperature reached in the fossil 
HII region models does not depend on which rate coefficient 
is used for reaction 13.) 

The effect of implementing the new rate for associa- 
tive detachment is less dramatic; note, however, that the 
uncertainty associated with this rate coeffici ent has been re - 
duced thanks to the recent study by iKreckel et al.l l|2010h . 
Using the systematic uncertainty given by the authors, we 
implement the thermal rate at the ± la levels, and find the 
following ranges for the critical flux. In the fossil HII region: 



accordingly, we restrict the discussion here to focus on two 
examples of thermal rate coefficients that have recently been 
updated, and how their revised values affect our results. 



7 I Glover et al. I j2006h have suggested that this rate is in fact erro- 
neously small, perhaps due to typographical errors in the source. 
Nonetheless, it is widely used in studies of star formation in 
metal— free gas. 
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Jcrit.HD = 3.85 ±0.15 x 10~ 3 [n/10 2 ]. In T vir > 10 4 K halos: 
Jcdt.HD = 3.7±0.2 x 10 _1 and 1.2±0.3 x 10~ 2 when illumi- 
nated by T5 and T4 spectral types respectively. For a thor- 
ough discussion of how this new rate coefficient compares to 
previous calculation s and its imp a ct on H2 chemistry, the 
reader is referred to lKreckel et al.l l|2010l 'l. 

Two final notes on the chemical network are on order: 
first, we have not included in this study the formation of H2 
via a three-body reaction, by which metal-free gas can be- 
come fully molecular at n > 10 s cm -3 , because our models 
do not include these high-density regimes. Lastly, as men- 
tioned above, the primordial value of D/H varies in the lit- 
erature by a factor of ~ 2, and this in itself may have con- 
sequences for the role of HD-cooling in metal-free gas. 
Fragmentation and Characteristic Protostellar Masses. Fi- 
nally, we emphasize that all our conclusions in this paper 
are based on the thermal history of a gas cloud, and how 
this is affected by a UV flux. In order to make realistic pre- 
dictions for the fragmentation, and the ultimate sizes of stars 
forming in a collapsing halo, fully 3-D hydrodynamical sim- 
ulations would be required. While robust HD-line cooling 
(or lack thereof) could have a notable impact on the char- 
acteristic stellar masses in the earliest dwarf galaxies, the 
process of fragmentation is yet to be fully understood, and 
thus physical ingredients such as the minimum gas temper- 
ature may not directly translate into the actual mass of the 
protostar that ultimately forms (see, e.g. JClark et al.ll2010l 
for recent results on fragmentation in metal-free gas, and in 
particular the importance of turbulence). 



5 SUMMARY 

We have demonstrated that HD-cooling in primordial 
gas can be suppressed by a relatively weak external LW 
background, with an intensity on the order of J21 ~ 
10 _3 (n/10 -2 cm -3 ) in const ant- density "fossil ionized" gas, 
and J21 ~ 10~ in shock- or photo-ionized gas collapsing 
into halos with virial temperatures greater than ~ 10 4 K. 
These critical intensities are lower than the expected mean 
UV background at z ~ 10 — 20, suggesting that HD-cooling 
is likely unimportant in most proto-galaxies forming near 
and just prior to the epoch of reionization. We conclude 
that an "HD-mode" of star formation was not as prevalent 
as previously thought. 

On a more technical note: we have also found that the 
negative feedback of the LW background is mediated via the 
abundance of molecular hydrogen, which is dissociated by 
the same radiation in its Lyman and Werner bands. Direct 
HD photo-dissociation is comparatively less important, al- 
though we find that in regimes of less effective self-shielding, 
it can regulate the minimum temperature of the gas. Finally, 
we have provided fitting formulae for the effects of HD and 
H2 self-shielding, shielding of both species by HI, and shield- 
ing of HD by H2, which we hope will be useful in other future 
studies. 
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